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Abstract 

It is shown that the two-body character of the interaction in a many-body system gives rise 
to specific correlations between the components of compound states, even if this interaction is 
completely random. Surprisingly, these correlations increase with the increase of the number 
of active (valence) particles. Statistical theory of transition amplitudes between compound 
states, which takes into account these correlation is developed and tested within the framework 
of the Two-Body Random Interaction Model. It is demonstrated that a feature, which can be 
called "correlation resonance" , appears in the distribution of the transition matrix amplitudes, 
since the correlations strongly reduce the transition amplitudes at the tails and increase them 
near the maximum of the distribution. 
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1 Introduction 



Measurements of parity non-conservation in neutron capture by Th nucleus gave a surprising 
result: in spite of a natural assumption of a random character of matrix elements of the weak 
interaction between compound states, the effect was found to be of the same sign for all observed 
resonances. Possibly, this means that the strongly fluctuating matrix elements of a weak pertur- 
bation between "chaotic" states of a compound nucleus (which was, in fact, the first example of 
quantum chaotic system) are essentially correlated. Thus, the problem of correlations between 
components of compound states is of great importance both for theory and applications. 

An attempt to study these correlations has been undertaken in p|, ^, It is natural to expect 
that some correlations may appear if the number of independent parameters in a Hamiltonian 
matrix is substantially smaller that the total number of the Hamiltonian matrix elements. The 
simplest example is given by the model of a random separable interaction p (see also @]), 

Hij = ei5ij + gViVj, i, j = 1,2, . . . , N, (1) 

where are the unperturbed energies and Vi are random variables distributed, e.g., according to 
a Gaussian law. As one can see, the number of independent parameters Vi is equal to A^, while 
the number of the Hamiltonian matrix elements Hik is A^^. It was shown that this model displays 
very strong (~ 100%) correlations between eigenvectors with close energies, despite the random 
character of the interaction: (viVj) = Sij{vf). 

Such correlations cannot appear in models described by full random matrices, like those of 
the Gaussian Orthogonal Ensemble (GOE). It was always obvious that such models possess some 
unphysical features, e.g., the semicircle level density, however, they seem to give a very accurate 
description of the correlations and fluctuations of energy levels (see, e.g., [Q). The important 
point is that in real many-body systems the basic interaction is a two-body one. This means 
that the number of independent parameters determining the n-particle Hamiltonian (two-body 
matrix elements) is much smaller than the number of the Hamiltonian matrix elements. Taking the 
two-body matrix elements as Gaussian random variables, a model called the Two-Body Random 
Ensemble (TBRE) was introduced in 0, § (see also references in 0). This model looks in principle 
much more realistic than the GOE. In particular, it has a Gaussian form of the level density, which 
is in good agreement with various nuclear shell-model calculations for "realistic" interactions in 
the finite basis (see, e.g., 0). The TBRE does not allow a deep analytical treatment, however, 
numerical modeling showed that its level fluctuation properties are very close to those of the GOE, 
although some differences were noticed Therefore, in some respect, the two-body nature of the 
particle interaction does not reveal itself in the level statistics. In other words, level fluctuations 
are insensitive to the details of the interaction between particles, provided the latter is large enough 
to cause strong mixing of the basis states. 

Another model which seems to be more physical then the GOE was proposed a while ago by 
Wigner to describe compound nuclei. It was suggested that the Hamiltonian matrix Hij has 
a banded structure, i.e., all matrix elements with \i — j\ > b are zeros {b is the bandwidth). The 
matrix elements inside the band are random numbers with zero mean and fixed variance, Hij = 0, 
Hfj = V'^, except those on the main diagonal, which monotonically increase. Ha = iD. In this 
model all eigenstates are localized in the unperturbed basis {V = 0). In the non-perturbative 
regime, 1 < V/D < \/b the strength function which describes the localization in the energy space 
(also called the local spectral density of states), has a characteristic Breit- Wigner form with a width 
r = 2t:V'^ /D within the band. This shape is in agreement with nuclear data |TT|, and with the 



calculated localization properties of chaotic eigenstates in the rare-earth atom of Ce |12[. It has also 
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been shown in that the Hamiltonian matrix which produces the dense spectrum of atomic 

excited states in Ce, is sparse and has a band-hke structure, although the edges of the band are 
diffuse. We should also mention that a number of analytical results on the localization properties 



of such Band Random Matrices (BRM) has been recently obtained in |jT^, |T5[ (see also [|T6[ and 
references therein). However, the off-diagonal Hamiltonian matrix elements in this BRM model are 
independent random variables, thus this model is void of any possible correlations related to the 
two-body interaction between particles. 

In this paper we show that there are quantities (transition amplitudes, or transition strengths) 
for which the underlying two-body interaction is of crucial importance. We show that the such 
interaction gives rise to specific correlations between the components of eigenstates, which are very 
essential for the distribution of transition strengths. Our results are obtained in the framework 
of the Two-Body Random Interaction Model (TBRIM) recently proposed in |]T^ for the study 
of various physical problems related to such complex many-body systems as heavy atoms, nuclei, 
metallic clusters, etc., which display quantum chaotic behaviour. Being in some aspects similar 
to the TBRE, the TBRIM is simpler in the sense that it abandons all restrictions imposed by the 
conservation of the angular momentum which makes it closer to the embedded GOE [|18|. On the 
other hand, the non-degenerate spectrum of the single-particle orbitals the TBRIM is based upon, 
generates a realistic level density and leads to a band-like structure of the Hamiltonian matrix. 

We should mention that there was quite a number of earlier works where strength distributions 
were studied using statistical spectroscopy methods and nuclear shell- model calculations (|l^, see 
also review 0). These methods are based on the calculation of distribution moments, which are 
given by traces of products of the operators in question and powers of the Hamiltonian over the 
model finite-dimensional space of the problem. Since the calculation of traces does not require 
knowledge of eigenstates, the question of correlations within eigenstates which is of prime impor- 
tance for the present work, has not been addressed in those studies. We must add that statistical 
spectroscopy methods emphasize and employ a particular importance of Gaussian spreading of 
many-particle configurations, and features like Breit-Wigner localization either do not appear, or 
are neglected (together with the interaction between configurations) in that formalism. All in 
all, it is unfortunately very difficult for the present approach to make contact with those results. 
Comparing the two approaches we should say that at first sight ours does not look as rigorous 
and mathematically advanced as the other one, as it appeals to some heuristic arguments and 
uses rather simple mathematics, e.g., perturbation theory. However, we believe that supported 
by numerical experiments, our method can give a deeper insight and a more physical picture of 
transitions between and correlations within the chaotic eigenstates in complex many-body systems. 

In Section 2 of this paper we show how the basic two-body interaction results in the correlations 
between the Hamiltonian matrix elements, eigenstate components and transition amplitudes. In 
Sec. 3 we check, whether the effects found in Sec. 2 could lead to some correlations between 
transition amplitudes coupling different pairs of states. Section 4 presents a brief outline of a 
statistical approach to the calculation of transition strengths; the analytical results are checked 
there against numerical ones obtained in the TBRIM. Finally, in Sec. 5 we study the spreading 
widths of the many-particle basis states. 
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2 Correlations between eigenvector components induced 
by two-body interaction 

Let us consider the basic ideas of the TBRIM. In this model, n Fermi-particles are distributed 
among m non-degenerate orbitals. In doing numerical experiments, we assume, as in fl^, that the 
energies of the orbitals are given by the simple expression: 

^a = do(^ + ^ , a = l,2, ...,m. (2) 

However, the analytical treatment presented below does not depend on a particular form of ■ 
Many-particle basis states \i) are constructed by specifying the n occupied orbitals. The energy Ei 
of the basis state equals the sum of the single-particle energies over the occupied orbitals. The total 
number of the many-particle states in the model is = n\{m-ny. ~ exp(?7, \n^ + {m—n) In The 
latter estimate relates to large m and n and shows that iV is exponentially large for n, m — n ^ 1 . 

The number of independent parameters of the many-body Hamiltonian is given by the number 
of different two-body interaction matrix elements V^/j-^^ and equals = m^(m — 1)^/2 . Due to the 
two-body character of the interaction, the Hamiltonian matrix element Hi^ = {i\ H \j) is non-zero 
only when \i) and |j) differ by no more than two occupied single-particle orbitals. As a result, the 
number K of the non-zero matrix elements Hij is given by 

K = N{Ko + + K2) , 

Kq = 1 , Ki = n{m — n), K2 = -n{n — l){m — n){m — n — 1) , (3) 

where Kq, Ki and K2 are the numbers of the Hamiltonian matrix elements coupling a particular 
basis state i to another one, j, which differs form i by the positions of none, one and two particles, 
respectively. Therefore, for n,m — n ^ 1 we have N2 K <^ N'^, i.e., the Hamiltonian matrix is 
essentially sparse and, in a sense, strongly correlated. 

To see the correlation between non-zero matrix elements, let us consider a pair of basis states 
\i) and |j) which differ by the states of two particles, for example, the state \j) can be obtained 
from the state \i) by transferring the particles from the orbitals a, (3 into the orbitals 7,5. For 
all such pairs, the Hamiltonian matrix elements are the same, Hij = Va/Bjs (or, strictly speaking, 
Hij = ^Va/B'ys, due to Fermi statistics). It is easy to calculate the total number Neq of the matrix 
elements Hij equal to Vajj-ys, using the fact that the remaining n — 2 particles can be arbitrary 
distributed over m — 4 orbitals, 

N - (4) 

(n-2)!(m-n-2)! • ^ ^ 

For basis states \i) and |j) which differ by the state of one particle {a — (3) the matrix element 
Hij equals the sum of the n — 1 two-body interaction matrix elements, Hij = J2'y Va-yp-y (the index 
7 runs over the rest n — 1 occupied orbitals). In this case Hij for different \i) and \j) (with fixed a 
and /5) do not coincide, but may contain identical terms Va-yp^^ , i.e., they are also correlated. 

The eigenstates \ni) of the model are determined by their components C\^^^ with respect to the 
many-particle basis states 

\ni) = Y.C?'^\^ , (5) 
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and can be found by solving the Schrodinger equation, 

= E("^)cf^^ . (6) 

j 

If the perturbation V is strong enough, the exact eigenstates |?t,i) are superpositions of a large 
number of basis states. As is known, strong mixing of basis states in the exact eigenstates (com- 
pound states) occurs locally within some energy range, Ei — E^"'^^ < T, where T is known as the 
"spreading width". It can be estimated as F N^D, where D is the local mean level spacing for 
many-particle states and Nyj is effective number of basis states represented in a compound state. 
This number is also known as the number of "principal components" . These components give the 

main contribution to the normalization condition C'f"^'' = 1 for the eigenstate |ni). Formally, 

we can estimate A^^ as the reciprocal of the inverse participation ratio, ^ . 

It is rather straightforward to show that the correlations between Hij result in correlations 
between the components C\^^\ Indeed, let us multiply the Schrodinger equation by the coefficient 
Ci"^) and sum over rii. Using the orthogonahty condition C*i"^^C*j"^^ = one obtains 

if,. = ^cf^)E('^Oc("i) . (7) 
"1 

In what follows, we assume that the matrix elements of the two-body interaction V are random 

variables with the zero mean, therefore, H^j = for i ^ j. In this case one can get Ct'^Ct^ = 
where the line stands for averaging over different realizations of V. However, if matrix elements of 
the Hamiltonian are correlated, HijHki ^ 0, the components of different eigenvectors |ni) and |n2) 
are also correlated, since 



nin2 



The latter relation shows that '^Cj"^^Ci"^^C/"'^ ^ 0. 

The above conclusion has important consequences. Let us consider a single-particle operator 

^ = E «l«/3^«/3 = E P"/3 Map (9) 
a/3 af} 

where aj^, and are the creation and annihilation operators. It is convenient to express the matrix 
elements of M in terms of matrix elements of the density matrix operator p^/j = a^^a^ which 
transfers a particle from the orbital (3 to the orbital a. One can see that the matrix element of M 
between compound states, 

(rill M 1^2) = E^"/3 (^i|P"/3 1^2) = 

a/3 

= E^a/3E^i"^^^|p«/3b)cf^ (10) 

a/3 ij 



has the zero mean due to the statistical properties of the components, i.e., pap \n2) = 0. Since 
the summation over the orbitals a,f3 in Eq.(^) is independent from the averaging over different 
realizations of V, in what follows we consider the simplest case of M = pap- The variance of the 
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matrix element of pa/3 between the two compound states is equal to 



{nil pap\n2) {n2 \ p^a\ni) 



ijkl 

ri(nin2) , c(nin2) 
O., "I- o„ , 



where we separated the diagonal and non-diagonal contributions to the sum (|np, 



S. 



(nin2) 



ik 



K) 



a 



'a/3 



q{nin2) 



i^j,k^i 



(12) 



(13) 



Note that the diagonal term S"^"^"^^ is essentially positive and can be easily estimated (see |T2|, ^ 



and Sec. 4 below), while the non-diagonal term S'^"'^"^) is of our main interest. If the eigenstates are 
completely "random" (different components both inside each eigenstate and of different eigenstates 
are uncorrelated) , the correlation sum Sc is equal to zero and the variance is determined by the 
"diagonal" sum 5*^ (this assumption have been used in the previous calculations of matrix elements 
between compound states in P, |l^, ^ 0)- However, we show below that in many-body system 
these two terms are of the same order, Sc ~ 5"^, even for the random two-body interaction V . 

The TBRIM allows one to investigate various properties of chaotic many-body systems taking 
into account the two-body nature of the interaction between particles. In the previous papers 



12, r7| there were indications that the diagonal approximation is not completely accurate for the 



computation of the variance of matrix elements of perturbation. In order to study this effect in 
detail, we have performed numerical experiments with TBRIM for the parameters corresponding 
to the model calculations of the Ce atom [|1^, 0. We take the number of particles n = 4, the 
number of orbitals m = 11, the spectrum of the single-particle orbitals is determined by do = 1, 
and the Gaussian random two-body interaction is given by Vv^ = 0.12. As a result, the size of the 
Hamiltonian matrix Hij is N = 330. The calculation of the matrix elements between compound 
states in this model gave a remarkable result. In Fig. 1 we present the "experimental" value of 
(see Eq. (|rT|)) together with the diagonal contribution (|T2|) (see Fig. la), and the ratio R = Sc/ Sd 
in Fig. lb. Figure la reveals a systematic difference between the diagonal approximation and exact 
expression ([TT|), and Fig. lb shows that non-diagonal term Sc is of the same order as Sd which 
clearly indicates the presence of correlations. 

Below, we show how these correlations emerge in the non-diagonal term Sc- First, note that for 
a given i the sum over k in Eq. (|1^ for Sd contains only one term, for which \k) = a^i^aa\i) = 
determined by transferring one particle from the orbital a to the orbital (3 in the state \i) (hereafter 
we will use the notation i' to denote such states). Accordingly, the index i runs over those states in 
which a is occupied and f3 is vacant. For such i and i' the matrix element (^l pap = 1, otherwise, 
it is zero. Therefore, in fact, the sum in ( P^ is a single sum, with a number of items less than A^, 



s 



(nin2) 



{ni) 



c: 



("2) 



(14) 



where the sum X)/ runs over the specified i. Analogously, Eq. (|T^ can be written as the double 
sum over i and j specified as above. 



q{nin2) 



E" r^{ni) ^{rii) ^{n2) ^{n2) 
L^j Oj/ Lyji , 



(15) 
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where j' is a function of j, = ajjOalj). Note that the energies of the basis states and their 
primed partners are connected as Eii — = ep — ea = Eji — Ej. 

One can expect that maximal values of the sum (JT^) and, possibly, ([T^) are achieved when 
C"s are principal components of the eigenstates. This means that mean square of the matrix 
element |(ni| pa/3 1^2) |^ is maximal when the operator pajj couples the principal components of the 
state \ni) with those of |n2), i.e. for ^'("i) — _E("2) ~ = _ p^j- from the maximum 

> r) a principal component of one state, say, rii, is coupled to a small 



component k of the other state n2 



E. - 



> r). The latter case is simpler to consider 



analytically, since the admixture of a small component in the eigenstate can be found by means of 
perturbation theory. This approach reveals the origin of the correlations in the sum Sc, Eq. (plSf). 
For example, if Cj"^'' is a small component of the eigenstate rii, then it can be expressed as a 
perturbation theory admixture to the principle components. If Cf"^^ is one of the latter, then there 



is a term in the sum ([T5|) , which is proportional to the principal component squared. 
Indeed, there are three possibilities: 
(i) ^) and Cf'^ 
small components. Then, one can write 



(ni) 



and Cp> are among the principal components, and C. 



(ni) 



and Cjl"^^ correspond to the 



("i) 



- E, 



E 



Hjp 



E P 



(16) 



c: 



("2) 



{t'\H\h2) 



E 



Hi'q 



n{n2) 



(17) 



Tilde above the sums indicates that the summations runs over the principal components only. 
The "coherent" contribution to the sum Sc in Eq. ([15|) is obtained by separating the squared 
contributions of the principal components in the sums in S'^"'^"^'* (i.e., p = i, q = j') 



Hi'j'Hji 



(E("2) -Ei'){Ei'^^) - Ej) 



a 



("i) 



("2) 



:i8) 



Taking into account that for the principal components we have E^ «i ^^('^i) and Ej/ ^ E^"'^\ we 



can replace the energies, Eii 
contribution to 



^("■1) _)_ CO and Ej £'("2) — upa, and thus obtain the following 



(^£;(".2) _ _E'(rii) 



(1) 



(2) 



Hi'j'Hij 



(19) 



(ii) C,-"^^ and Cf"^^ correspond to the principal components, C^-"^^ and C-^^^ correspond to the 



small components. Then, the result is the same as (|T^. 



iii) Cf''^'' and C-'^^'' are principal components, C^?^^ and (7-7^'' are small components (or, C^, 



(n2) 



and d?'^ 



y are principal components, Cj-"^'' and Cj"^'' are small components). In these cases there 



are no coherent terms in the sum for Sc in Eq. (0). This follows from the fact that for chaotic 
eigenstates the mixing among the principal components is practically complete, which makes them 
to a good accuracy statistically independent. 

> r, one obtains 



Thus, far from the maximum. 



^("2) _ ^(ni) _ 

2 



(^_£;(n2) _ _£;(rii) 



x2E 

'^t3a) ij' 



(1) 



(2) 



Hi'j'Hij 



(20) 
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A similar calculation of the diagonal sum S^J'^^^\ Eq. (|T2|), yields 



S. 



(nin2) 



X 



("i) 



("2) 



Hh + EE 



(ni) 



(«2) 



(21) 



The two terms in square brackets result from the contribution of principal i and small i' components 
in Eq. (0), and vice versa. From Eq. pO| ) we see that = if //j'j'ifij- = 0. However, there 

is nearly a hundred percent correlation between these matrix elements. Indeed, the basis state i' 
differs from i by the location of only one particle (the transition from the orbital a to (3), and the 
same is true for j' and j. 

Let us estimate the relative magnitudes of Sd and Sc- First, consider the case when \i) and \j) 



i). In this case Hij = 



UlfllU2tJ.2 ■ 



V,. 



Since the basis states 
Hij (note that 



differ by two orbit als, |j) 

and must differ by the same two orbitals, we have Hiiji 
^1, fi'2 7^ C(,P, since both states \i) and \j) contain a and do not contain /3, whereas and 

contain j3 and do not contain a). Therefore, the averages over the non-zero matrix elements 
between such pairs of states are: HijHi'j/ = Hfj = Hf,-, = V"^. 

Now, let us consider the case when \i) and |j) differ by one orbital |j) = al^ai,^\i). In this case 
the Hamiltonian matrix elements are sums of the n — 1 two-body matrix elements, 

n-2 

Hij ^ ^ ^uifj,U2fJ, ~l~ '^V\au2a ; 
n-2 

^i'j' = E ^i'iMi'2M + ^i^i/3i'2/3 • 

The sums of n — 2 terms in Hij and Hi/j/ coincide, the difference is due to the one term only (orbital 
a is replaced by the orbital P). Thus, 



HijHi'j' 



{HijY = {H,,>Y = {n- 1)V' 



where we took into account that K. 



The contributions of one-particle and two-particle transitions in Eqs. ( PP]) and (^Tj) representing 
Sc and 5*^ respectively, will be determined by the numbers of such transitions allowed by the 
corresponding sums. For the single-prime sums in Eq. (^) these numbers are proportional to Ki 
and K2, Eq. (^. In the double-prime sum in Eq. ( pOf ) these numbers are proportional to Ki 
and K2, the numbers of the two-body and one-body transitions i — >■ j', in the situation when one 
particle and the two orbitals {a and (3) do not participate in the transitions. These numbers can be 
obtained from Eq. (^) if we replace n by n — 1, and m by m — 2, so that Ki = (n — l)(m — ri — 1), 



[n- 



l)(n-2)(m-n-l)(m-r2-2)/4. Finally we obtain that at 



> r the 



contribution of the correlation term to the variance of the matrix elements of pa/s can be estimated 
in the ratio as 



R = 



Sc 



(n 



2)K^ + K2 



{n - l)K^ + K2 
{n — 2){m — n — l){m — n + 2) 
n{m — n){m — + 3) 



(22) 
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This equation shows that for n = 2 we have Sc = 0, which is easy to check directly, since HiijiHij = 
in this case. For n > 2 the correlation contribution Sc is negative at the tails of the strength 
distribution. This means that it indeed suppresses the transition amplitudes off-resonance (see 
Fig. 1). For n,m — n ^ 1 the ratio R is approaching its limit value —1. It is easy to obtain from 
Eq. (ID that for m - n > 1 

^—^ = l + R^- -. 23 

bd n[m — n) 

Thus, surprisingly, the role of the correlation contribution increases with the number of particles. 

For the numerical example shown in Fig. 1, n = 4, m = 11, one obtains: R = —0.39, which 
means that the correlation contribution reduces the magnitude of the squared matrix elements 



> r ). The ratio found 



between compound states almost by a factor of 2 (for 
numerically is i? ~ —0.45 (Fig. lb, n2 =150-250; larger n2 are probably too close to the boundary 
of the matrix for R to remain constant). 

We would like to stress that the role of the correlation term does not decrease with the increase 
of the numbers of particles and orbitals. This prediction is supported by Fig. 2, which shows 
the behaviour of the squared matrix element and its diagonal and correlation parts for n = 7 and 
m = 14, = 3432. One can see that the suppression of the matrix elements due to the 
correlation term at the tails is even stronger than that in Fig. 1 [the numerically found ratio is 
R ^ —0.7 Ys R = —0.55 obtained from Eq. (p2D]. The correlation contribution should be even 
more important in compound nuclei, where N ~ 10^ . This case can be modeled by the parameters 
n = 10, m = 20; then we have R = —0.66, or, equivalently, {Sd + Sc)/Sd = 0.34, which means that 
the correlations suppress the squared element between compound states by a factor of 3 (far 
from its maximum). 

It is worth emphasizing that the existence of correlations due to the perturbation theory admix- 
tures of small components to the chaotic eigenstates, which leads to a non-zero value of Sc ([I3|), is 
indeed non-trivial. For example, if one examines the summand of Eq. (|15]) as a function of i and 
j, it would be hard to guess that the sum itself is essentially non-zero, since positive and negative 
values of C^^^^ Cj^^^ C^J^'^'' C^?^^ seem to be equally frequent, and have roughly the same magnitude, 
see Fig. 3. 

Since Em '5'^'''"'^ = E„2 5'^"^"") = (see below), the suppression of JP at the tails should 
be accompanied by correlational enhancement of the matrix elements near the maximum (at 



— uJi3a\ < r). Thus, we come to the important conclusion: even for a random two-body 
interaction, the correlations produce some sort of a "correlation resonance" in the distribution of 
the squared matrix elements M^. One should note that this increase of the correlation effects in 
the matrix elements of a perturbation can be explained by the increased correlations between the 
Hamiltonian matrix elements when the number of particles and orbitals increases {N/n oc e"). 
Now we can estimate the size of the correlation contribution Sc near the maximum of the the 



distribution (at 



< r). First, we show that after summation over one of 



the compound states, the correlation contribution vanishes. Indeed, 

"2 "2 ijtj.k^l 

where we take into account that the sum over n2 in the expression above is zero for k ^ I. Therefore, 



the negative value of Sj^^'^"'^^ at 



> F must be compensated by its the positive 
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value near the maximum. The sum rule (|2^) allows one to make a rough estimate of Sc near the 
maximum of Sd (and M^). 

Let us assume that = RmSd at lE^"^) _ ^("O _ ujpj\ < V/2, whereas = RtSd at [E^^^) _ 
- C^/3a| > r/2 [Rt is given by Eq. (^)]. The distribution of S^^^^^ can be reasonably 
approximated by the Breit-Wigner shape (see Sees. 1 and 4), 



Q(ni"2) 



A 



+ rV4 ' 

where E = E^^'^^ — i^^^i) — u;^^, and F = + F^^- The sum rule (|^) implies that 



(25) 



Rr, 



r dE r 



r/2 



E2 + r74 



(26) 



Since the two integrals in the above equation are equal, we have Rm = —Rt- Thus, near the 
maximum the correlation contribution Sc is positive and enhances the squared matrix element with 
respect to the diagonal contribution. 



Sd + Sc 

~^7~ 



l + R„ 



m 



n{m — n) 



(27) 



Comparing the values of the ratio Sc/ Sd at the maximum and at the tail in Fig. lb (n = 4, m = 11), 
one can see that indeed, Rm ~ —Rt- For larger n and m the correlation enhancement factor 
asymptotically reaches its maximal value of 2. The numerical example in Fig. 2 (n = 7, m = 14) 
shows the enhancement of with respect to Sd at the maximum even greater in size than that 
predicted by Eq. (|27|) . This is not too surprising since in Eqs. (p5D-(p7D we estimated the average 
value of Rm over an interval AE ~ F around the maximum rather than the peak value at the 
maximum. 

A similar estimate of Sc near maximum can be obtained by the direct calculation of the small 
component contribution to Sc (|1^)- On an assumption that there are no correlations between 
principal components of compound states we can separate the contribution of small components. 
For example, in the resonance situation, i?^"^) _ ^(m) ^ u/^a, if the components Sj"'^'' and Sjr^^ 
are small {\Ej — E^"-'^'^\ > F, and consequently, \Eji — E^'^^^\ > F), then they contain contributions 



proportional to the principal components and C^v [see Eqs. (|T6D, (|17D]. Analogously, S. 

and may be among the principal components, and then the small components Cf"^^ and Cf"^^ 
will contain correlated contributions. Thus, we have the following estimate. 



(^i) 



small j 



C! 



("2) 



HijHi'j' 



(^(ni) _ Ej){E^'^^) - Ej 



(28) 



Since Ej' 



Ej = Ei' — Ei c:i _e("2) _ ]^{ni) ^Qj, ^YiQ principal components i and z', (E^^^^ — Ej) 
and (£^("2) _ ]^.,'^ ^]-^g denominator have always the same sign, and Sc is positive (recall that 
HijHi'j' > 0). Expression ( pSD can be estimated using the well known formula for the spreading 
width, F = 27THfj/D, where D is the mean level spacing for the many-body states. This yields 



Sc ~ Sd, in agreement with the previous estimate (P?]). 
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3 Correlations between transition 
amplitudes 

We have shown that correlations between eigenvector components in a system with a two-body 
interaction between particles must be taken into account when calculating the variance of a matrix 
element between compound states M^. Another question is whether the above correlations between 
eigenstate components lead to correlations between different matrix elements, 

M„j„2M„2„3 = {ni\pap 1^2) (^2! P/3a \n^) 

= E Ct'^C^'^'^Ct^Ct'^ {^\ Pa, \j) {k\ \l) . (29) 



Our analysis shows that the correlations of the type (|29| ) are absent, i.e., M^^njikf^^^g 



M„j„2 M„2„3 = 0). The result could be different if the principal components of different eigen- 
vectors (l^^i) and 1^3)) were correlated. This effect takes place in the separable interaction model 
P, H, but we have not found such correlations in the TBRIM. 

The absence of correlations between different amplitudes is confirmed by direct numerical ex- 
periments. First, we have studied the probability density of the matrix elements M„^„2 for different 
ni,n2 obtained for a number of realizations of the two-body matrix elements V^^^^. Since the 
variance of M„j„2 depends on rii and n2, the probability density of M„j„2 has been obtained by 
normalizing each matrix element M„^„2 its root-mean-squared value which was calculated by 

averaging over the realizations of Vafs^s- The resulting probability density V ^M„j„2/ \J ^nm-i^ ; 

averaged over A^^ = 5 realizations of Vq,/3^5, turns out to be quite close to Gaussian (Fig. 4). This 
result follows from the fact that each matrix element between compound states is the sum of a large 
number of random (or, almost random) terms, see Eq. (p!OD, so that the Central Limit Theorem 
applies. Therefore, the correlations found in the previous section, do not show up, unless more 
complicated correlations involving different components of the same eigenstate, like those in Eqs. 
(H) or ([TT|), are probed. 

To check whether some correlations between different matrix elements (PP| ) exist, we have plotted 
the matrix element {ni\pa(3\n2) versus another one, (ni|po/3|n2 — 1), where |n2 — 1) is the eigenstate 
immediately preceding |?7,2), for some fixed ni, n2, a, and (3, obtained from = 387 different 
Hamiltonian matrices (Fig. 5a). The latter were generated by using different random realizations 
of Vafs^s- Detailed analysis of the distribution of the points in this figure does not reveal any sort 
of correlations. 

The next question is the existence of correlations between matrix elements M„j„2 ^'^d lVnin2 of 
different operators for the same compound states |?t,i) and \n2). If the expansions of these matrix 
elements [see Eq. (|10|)] contain identical matrix elements of the density matrix operator p^fs, such 
correlations, in principal, do exist: 



Mn.n^Wn^n, = E ^«/3^/3a (^1 1 P"/? I^s) ' 7^ . (30) 

A more complicated question is whether the matrix elements of different "elementary" transition 
operators pais and p^s are indeed uncorrelated [as we assumed writing Eq. (|30D]. The product of 
such two matrix elements can be presented in the form 

{nMn2){n2\ps,\n,) = Ecf"^^C7f^^Cf^)C^)(^|p«^|A:)(/|p,,|j) (31) 



ijkl 



c 



(ni) 



c. 



("2) 



2 



'^p {E^^^) - E,){Ei^^) - E,) ' ^^^^ 
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where the last expression is written for rii and n2 far from the maximum _ ^(^i) _ ^^^^ > 

r, \E^'^^'' — iJ^^i) — LJsyl > r), and = ajjOali), = aja^|j). It can be shown that in our model 
HijHiijii oc ^a-y^ps- Therefore, there are no terms in the expression (|32D which would give non-zero 
contributions, and the average of (|3l|) is zero. The absence of correlations in this case is illustrated 
by Fig. 5b, where numerical data obtained in the TBRIM are presented. As in the case of the 
matrix elements between different pairs of compound states, no correlations can be seen between 
the matrix elements of different transition operators. 



4 Statistical description of the transition amplitudes 

In this Section we use the TBRIM to test the validity of the statistical approach to the calculation 
of transition amplitudes between compound states of complex systems developed in [|, . In 

what follows we first outline the main ideas of the statistical approach. The variance of the matrix 
elements of an operator M (|^) between the compound states and 1^2) can be presented in the 
following form [compare with Eq. (PPP], 



IM, 



nin2 I 



a/3 



(33) 



where we have taken into account the result of the previous section that the average of the correlator 
(|3TD is zero unless 7 
to the calculation of 



a, 6 



(or M„^„2pV„2„,; 



is reduced 



p. Therefore, the calculation of jM^^n^l 

2 can be presented as the sum of the diagonal 



It was suggested in Sec. 2 that \{ni\pai3\n2 
and correlational sums Sd and Sc, Eqs. (|T^)-([T^). Since we have already estimated the ratio 
{Sd + Sc)/Sd, it is enough to calculate only Sd, Eq. (|12D. Following let us replace the squared 



components 



C 



(«i) 



and 



("2) 



by their average values. 



a 



(ni) 



W 



C 



(n2) 



W 



(34) 



where the averaging goes as usual either over a number of realizations of the two-body interaction 
matrix elements (ensemble average), or over a number of neighbouring eigenstates (physical energy 
average); in the spirit of ergodicity the results are presumably the same. The function w is pro- 
portional to the strength function introduced by Wigner which is also called the local spectral 
density of states. Note that definition (|34D also implies that the mean-square contribution of the 
component i in the eigenstate rii is determined by their energies, Ei and (in fact, by their 

difference \Ei — E^"''^^\). For states localized in the given basis, w is a bell-shaped function with a 
typical width determined by the spreading width F. There is some theoretical and experimental 
evidence that it can be approximated by the Breit-Wigner formula, although its tails decrease faster 
then \Ei - (see references in Sec. 1). 

The diagonal sum now takes the form 



r<{"in2) 



ik 



w{Ei,E^''^^)w{Ek,E^^''>){iMk){k\pp^\z) 



(35) 



The summation over k for a fixed i includes only one state, \k) = Pf3a\i), with E^ = E^ + ojpa- On 
the other hand, we can write 



^{i\pci3\k){k\ppa\i) = {i\pappf}a\i) = {i\na{l - nf3)\i) 

k 



(36) 
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where and = a^^ap are the occupation number operators. Thus, we obtain, 

i 

The matrix element (i|na(l — np)\i) is equal to 1 if the orbital a is occupied and (3 is vacant in 
the basis state otherwise, it is zero. We used this fact earlier [Eq. (|T^) ] to reduce the summation 
to these states only. Now we proceed in a different way. Both w's in (^) are smooth functions 
of energy normalized as Y^i^iEi, E^"^^^) = 1. This allows one to replace the matrix element of 
no,(l — fip) by its expectation value. 



= E \ct'^\\^\na{^ - ^/3)K) - (^a(l - np))n, . (38) 

i 

The sign ~ above is a reminder that the left-hand side is the local average over the states |ni). 
Practically, when the number of components is large, the fluctuations of (na(l — n/3))„i are expected 
to be small. Now we can re- write Eq. (pTf ) in a form similar to Eq. (|T^), but without any restrictions 
on the summation variable i, 



It was shown in [|T^ that under some reasonable assumptions about the functions w one can 



introduce a "spread 5-function" (5(A), 

5(A) = D,' E ^(^- E(^^^)w{E., + cop^, 



J -^wiE,, E^-^^)w{E, + c.^,, E("^)) , (40) 



where A = E("2)-^("i)_c^^„, Di and D2 are local mean level spacings for the rii and n2 eigenstates. 
The function (5(A) is symmetric, its characteristic width is determined by the spreading widths of 
the eigenstates rii and n2, F ~ Fi + F2, and it is normalized to unity, / (5(A)(iA = 1, just as 
the standard (5-function. If w's have Breit-Wigner shapes, 5 is also a Breit-Wigner function with 
F = Fi -|- F2. The fact that SjP^"'^^ is proportional to the function (5(A) is a particular manifestation 
of the energy conservation for transitions between the quasistationary basis states ||2^ [if F — >■ 0, 
then (5(A) (5(A)] . Using Eqs. (0), (p9D, and ( ^0]) we can finally present the diagonal contribution 
to the variance of the matrix element M„j„2 in the following form 



l^n.n.l'diag = E l^apf R(l " ^/3))ni/^25(^^"^^ - E^^^^ - COp^) . (41) 
a/3 

This expression is apparently asymmetric with respect to the states rii and n2. By performing the 
calculation in a different way we can obtain 

^("^"^) = {np{l-n^))^^J2^iEk-upa,E^'^^^)w{Ek,E^-^^) , (42) 

k 

instead of Eq. (P^), and thereby arrive at a different formula for the variance. 



a/3 
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where the occupancies factor is now calculated for the state n2 (it represents the probability to 
find the orbital (3 occupied, and a empty). If the suppositions made in the above derivations are 
correct, the two formulae (^Tj) and should give identical results. 

In the present work we use the TBRIM to check the accuracy of the statistical approach de- 
scribed above. Figure Fig. 6a presents a comparison between the values of S*^"^"^^ as given by Eqs. 



(pOD, (|4^), and those from the initial expression (|T2D. Clearly, there is a good agreement between 
the three formulae. 



It is quite important for applications of the statistical approach (see 0, ^) that further sim- 
plifications be made by replacing the correlated occupancies product {hahp),^^ in Eq. (|38|) by the 
product of the two mean values, {na)nx{ni3)nx- This is definitely a valid operation when the num- 
bers of excited particles and active orbitals are large, so that the occupation numbers for different 
orbitals become statistically independent. Then one would be able to use the following relation 

(na(l - fip)) ~ n{t^)[l - n{ep)\ , (44) 

where n{ea) and n{ea) are the occupation numbers. (They can be calculated, e.g., using the 
Fermi-Dirac formula with an effective temperature, see see also |^ where the relation 

between thermalization and chaos is studied in nuclear shell-model calculations). The result of 
such simplification is shown in Fig. 6b, where the diagonal contribution (O) is again compared 



with the values obtained from Eqs. (p9D, (^21) , using approximation (]4^) . In spite of the fact that 
the TBRIM calculation included ra = 4 particles only, the agreement remains quite reasonable, the 
error being about 10%. To examine the quality of the approximation at the tails of the distribution. 



Fig. 7 shows the ratio of Sd as given by Eqs. (|39|), (^21) to the directly calculated diagonal term 
S'd, Eq. ([T2|) . The difference between Fig. 7a and Fig. 7b highlights the inaccuracy introduced by 
an additional approximation ( ^4|) for the occupation numbers. 

In order to make a more direct test of the validity of substitution (^^, we plotted in Fig. 8 the 
correlator (nQ,n/3)nj/[(nQ,)„j(n^)„J as function of ni. Consistent with the small number of particles, 
this correlator displays large fluctuations, however, its average value of about 0.8 is still rather close 
to 1. 



5 Spreading widths for different basis components 

In Section 4 when considering the statistical approach to the calculation of the variance of matrix 
elements between compound states, it was assumed that the spreading widths F are the same for 
all basis components. However, this question is not trivial. As is discussed in the literature, the 
spreading widths of components corresponding to different number of excited particles could have 
significantly different values. For example, in p2[ it is argued that two-particle-one-hole states 
{2p-lh) can lead to correlations between values of parity nonconserving effects []T|, if the spreading 
width of the 2p-lh states is two orders of magnitudes smaller than that of Ip states. In such case 
one might expect that in our model the spreading width showed a rapid decrease as function of the 
number of excited particles in the basis component. 

To study this question in detail, we have performed additional tests. In Fig. 9 the root-mean- 
squared spreading width Tj for all basis states \j) is presented for n = 6, m = 12 (A^ = 924). Here 
we use the following definition: 

T', = {j\{H~{j\H\j)f\j) , (45) 
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which can also be presented in several equivalent forms, 



= i^^j,, - (HnY (46) 

= E^S (47) 

= M2 - (Ml)' , (48) 

the last one relating Tj directly to the moments Mp of the strength function pw{E,j) = 



J p^{E,j) EME = Y: f {E^^'^y . (49) 



Equations (^B|)-(^S]) can be obtained using closure, J2j \ = J2n = 1; where |j) and \n) 

are the basis state and the eigenstate respectively, {j\n) = Cj"''. We should note that the r.m.s. 
spreading width is different from that introduced intuitively in Sec. 4 as the characteristic width 
of the strength function. For example, if the strength function has a Breit-Wigner form, its second 
and higher moments are infinite. In the Wigner's BRM the r.m.s. spreading width is determined 
by the bandwidth b as Tj = a/26V^, whereas the "Breit-Wigner spreading width" is Fbw = '^ttv'^/D 
(see Introduction). However, in a more realistic situation the strength function drops rapidly, its 
second moment is finite, and the difference between the r.m.s. F and Fbw is not large. 

From Figure 9 one can see that apart from small natural fluctuations, the r.m.s. spreading 
width is the same for all components. To exclude a weak dependence of Tj on the energy of the 
basis state j (boundary effects seen as rises of Fj at small and large j), we calculate the mean value 
F and the r.m.s deviation ST of the spreading width by averaging over j =50-874. The results are: 
F ^ 2.22, 6r ^ 0.12. The latter value shows that the fluctuations of the width are very small. 



This result is in agreement with computations made for the Ce atom |]T2[; similar results have been 
recently obtained in the nuclear sci-shell model calculations [§, ^ . The fluctuations of the width 
are small due to the large number of "decay channels" for each basis component (each component 
is coupled to many others by random interaction) [Q. Formally, this can be obtained from Eq. 
(^71). For example, in the TBRIM one obtains 

p = (n - 1)V^K^ + V^K2 = K^2V^ , (50) 

which shows that for large number of "independent decay channels", K12 = {n — l)Ki + K2 [see 
Eq. (^] the statistics of Fj is given by the distribution with K12 ^ 1 degrees of freedom, 
resulting in the ~ 1/ K12 decrease of fluctuations. More accurately, the relative r.m.s. fluctuation 



of the squared width (^51) is given by (5(F^)/vF2 = y2/_ft'i2. For n = 6, m = 12 Eq. ( pOD yields 

vT^ ^ 2.41, and the relative fluctuation is 0.07. These values are close to the numerical ones 
quoted above (the discrepancy is mainly due to the difference between the mean width and its 
r.m.s. value, and the corresponding difference in the fluctuations of Tj and F^). 

To check the independence of Fj of the number of excited particles in the component j, we have 
calculated the average spreading width F(p) for basis states with a fixed number of excited particles 
p, p = 1, 2, n — 1. In the numerical experiment shown in Fig. 6 all for p =1-5 proved to be 
approximately the same, F(p) 2.2. 

The above consideration shows that the statistical approach does not provide any support for 
the dependence of the spreading width on the number of excited particles. This indicates that 
the argumentation in favour of a strong dependence, based on different decay phase volumes for 
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different numbers of excited particles, seems to be incorrect. In our opinion, the difference in the 
spreading widths could appear as a result of specific dynamical effects. For example, this could be 
an influence of levels in other potential wells which appear at higher nuclear deformation, or due 
to an interaction with collective motions, such as rotations and vibrations. 

6 Conclusions 

The calculation of the mean square matrix element of an operator between compound states of a 
many-body system has been considered. We have shown that the two-body nature of the interaction 
between particles manifests itself in the existence of correlations between the components of the 
"chaotic" compound eigenstates. These correlations taken together with the correlations between 
the many-particle Hamiltonian matrix elements result in a relatively large correlation contribution 
to the mean square matrix element. The correlations exist even if the two-body matrix elements 
are independent random variables, as in the TBRIM. Such correlations can be understood in terms 
of the perturbative mixing of the distant small components to the principal components of the 
eigenstates. If the Hamiltonian matrix elements are random variables the correlations of this type 
vanish. 

One of the most interesting feature of the correlations found in our work is that they do not 
decrease with the increase of the number of excited particles or active orbitals. Thus, they must 
be taken into account when calculating matrix elements of a weak interaction between compound 
states in nuclei. 

Another feature concerns the shape of the distribution of the density matrix operator near its 
maximum. As one can see from Figs. la,2a, the correlations create a sharp spikelike form of 
the distribution, instead of a smooth Gaussian or Breit-Wigner form. With such sharp peaks, 
the strength function for any particular operator M can have the so-called gross-structure, due to 
many single particle transition terms in the expression (0). Without these specific correlations, 
the strength function would be much smoother and the gross-structure would not be seen. It is 
also interesting to note that there are very large mesoscopic-type fluctuations in the distribution 
near the maximum, depending on a specific (random) realization of the two-body interaction V . 
This fact is also the consequence of strong correlations. 

Our study also demonstrated that the spreading widths of different basis components are ap- 
proximately constant and fluctuate very weakly. In particular, we have not found any dependence 
on the number of excited particles in the component. 

The statistical approach to the calculations of such matrix elements has been tested in the 
present work with the help of the TBRIM. The numerical results obtained in this work support 
the validity of the statistical approach. The TBRIM has also enabled us to check that the matrix 
elements of different transition operators between a pair of compound states are uncorrelated, as 
are the matrix elements of a given operator between different pairs of compound states. 
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Figure 1. 



FIGURE CAPTIONS 



(a) Mean square matrix element ( pTj) calculated in the TBRIM for n = 4 particles and m = 11 
orbitals, a = 4, /3 = 5, as a function of the eigenstate n2 for rii = 55. Averaging over A^^ = 100 
Hamiltonian matrices Hij for different realizations of the random two-body matrix elements has 
been made. Dots correspond to the sum Sd + Sc while the solid line represents the diagonal contri- 
bution 5*^ only [see (p!2D]. 

(b) Ratio R = Sd Sd of the correlation contribution to the diagonal contribution. 
Figure 2. 

Same as in Figure 1, for n = 7, m = 14, a = 7, /? = 8. The data obtained for a single Hamilto- 
nian matrix of the size A^ = 3432; rii = 575. Note the increased role of the correlation contribution 
Sc. 

Figure 3. 

The distribution of the items of the sum (0), = '^C]"^^Cf for m = 55, n2 = 66, 
obtained in the TBRIM for the same set of parameters as in Fig. 1, averaged over A^,. = 100 
realizations of V. Indices i and j in the figure run over those 84 components in which a is occupied 
and P is vacant. 

a. Positive values. 

b. Negative values (absolute values). 
Figure 4. 



Probability density of the normalized matrix elements x = {ni\pai3\n2) / (^|(^i|Pa/3|^2)| j in 
the TBRIM for the parameters of Fig. 1. The histogram is obtained for A^^ = 5 Hamiltonian 
matrices. Solid curve is the normalized Gaussian distribution. 

Figure 5. 

(a) The matrix element yap = {ni \ [^2) for ni = 55, n2 = 67, plotted vs. the matrix element 
Xai3 = {ni\ Pa/3 1^3) for ns = 66; a = 4, /3 = 5 and other TBRIM parameters as in Fig. 1. The 
number of points in the figure is A^^ = 387. No evidence of correlations between x and y is present. 

(b) The matrix element ya/s = Pa/s {^2) for a = 4, /3 = 6 vs. Xa-y = Pa-y \n2) with 7 = 5 
for 111 = 55 and n2 = 66. Again, there is no indication of correlations between x and y. Note the 
difference in the vertical and horizontal scale due to the fact that for given ni and n2 the energy 
difference — £^("2) jg approximately in-resonance for the transition between a = 4 and /5 = 5 
and off- resonance for a = 4 and 7 = 6. 

Figure 6. 
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(a) The diagonal contribution to the mean-square matrix element as obtained from Eqs. (pQ]), 
(^H) (solid lines) in comparison with the direct calculation of 5*^, Eq. (0) (circles). The TBRIM 
parameters are the same as in Fig. 1. 

(b) Same as (a), with the occupancy factors in Eqs. (l39|), (^2|) calculated by means of Eq. (^^ . 
Figure 7. 

(a) The ratio Ra of the approximation represented in Fig. 6a by the solid lines, to the value of 



(b) Same as in Fig. 7b for the data of Fig. 6b. 



Figure 8. 



The correlator Qai3 = {nani3)ni/[{fia) minis) m] vs the eigenstate number rii. The TBRIM pa- 
rameters are the same as in Fig. 1. The average value of the correlator is about 0.8, which means 
that the correlations between the occupancies of different orbitals are not very strong. 



Figure 9. 



The spreading width Tj calculated as the r.m.s. deviation from the center of the distribution 



of the components 
corresponding to n ■ 



("i) 



for each basis state j. The data are obtained for one matrix Hij 
6 particles and m = 12 orbitals. 
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